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For a finite XY chain and a finite two-dimensional Ising lattice, it is shown that the paramagnetic 
ground state is adiabatically transformed to the GHZ state in the ferromagnetic phase by slowly 
turning on the magnetic field. The fidelity between the GHZ state and an adiabatically evolved 
state shows a feature of the quantum phase transition. 
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I. INTRODUCTION 

A quantum computer is a quantum system, so it could 
simulate quantum dynamics directly and efficiently [1]. 
When a quantum system undergoes a quantum phase 
transition (QPT) 0, induced by the variation of external 
parameters or coupling strength, its ground state changes 
dramatically and maybe, too, its entanglement. Entan- 
gled states, showing quantum correlations between sub- 
systems, are not only valuable resources in quantum in- 
formation processing but also important for understand- 
ing quantum many-body systems. So much attention has 
been paid to a study of entanglement of quantum many- 
body systems in ground states or at thermal equilibrium 
in connection with QPTs [1, However, the simula- 
tion of QPTs and the generation of entangled states with 
quantum computers are less explored. 

In this paper, we address whether QPTs can be sim- 
ulated with adiabatic quantum computation (AQC) and 
how the entanglement changes when the QPT takes 
place. We present a way to generate an entangled state 
of a quantum system which undergoes a QPT during the 
quantum adiabatic evolution. As prototypes of QPTs, 
we consider a spin 1/2 XY chain and a two-dimensional 
spin 1/2 Ising lattice. It is shown that a product state 
in the paramagnetic phase is adiabatically transformed 
to a Greenberg-Horne-Zeilinger (GHZ) entangled state, 
in the ferromagnetic phase, and vice versa. We shows 
the fidelity between the GHZ state and an adiabatically 
evolved state could be a good indicator to QPTs. For 
a two-dimensional Ising model, a two-dimensional GHZ 
state is generated via AQC. 



II. ADIABATIC QUANTUM COMPUTATION 

Let us start with a brief introduction to AQC 
Quantum computation can be implemented by the con- 
trolled dynamics of quantum states governed by the 
Schrodinger equation with a time-dependent Hamilto- 
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and quantum measurements. In a quantum circuit 
model, the evolution of a quantum state is decomposed 
into a series of single-qubit and two-qubit gates, which 
can be implemented by applying external pulses and by 
controlling the interaction between two qubits. On the 
other hand, AQC relies on the adiabatic theorem, which 
states that an evolved quantum system will stay at its 
instantaneous eigenstate if the time-dependent Hamilto- 
nian changes very slowly. If the Hamiltonian H{t) in 
Eq. ID) changes slowly, then the initial state |5'(0)), taken 
as an eigenstate state |(/3„(0)) of an initial Hamiltonian 
7?(0) = i?o at t = 0, evolves to \^{T)) = |(/?„(T)), 
an eigenstate of a problem Hamiltonian H{T) = Hp 
at t = T. The run time T is inversely proportional to 
the square of the minimum energy gap during the evo- 
lution. It is convenient to introduce the dimensionless 
time s = t/T with < s < 1. There are many ways 
to connect Hq and Hp smoothly as a function of s, for 
example, simple linear or nonlinear interpolations [7]. A 
general interpolation is given by H(s) = f{s)HQ-\-g(s)Hp 
where two functions f{s) and g{s) satisfy the boundary 
conditions, /(O) = .9(1) = 1 and /(I) = .9(0) = 0. It is 
known that a proper interpolation could reduce the run 
time of AQC. 



III. ENTANGLEMENT AND QPTS OF THE XY 
CHAIN VIA AQC 

The Hamiltonian of a spin 1/2 XY chain in a transverse 
magnetic field is written as 
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where N is the total number of spins, A the trans- 
verse magnetic field, and 7 the parameter for the de- 
gree of anisotropy of spin-spin interaction. Here, the co- 
ordinates X and z are exchanged for convenience as in 
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Ref. 0. The periodic boundary condition, (Tat+i = cri, 
is assumed. For 7 = 1, it becomes the Ising model 

Hi = ~ Sill ^t^i+i - ^ Z]j ^i- Foi' 7 = 0, it is caUed 
the XX model. 

Let us recall the ground state of the Ising model which 
undergoes the QPT at Ac = 1 H. If A > 1, the Zee- 
man term in Eq. ([2]) is dominant and the ground state 
is given by the product of eigenstates of af , called the 
paramagnetic state 
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where |+) = ^ (|T) + In the other limit of A < 1, 
the spin-spin interaction in Eq. ^ is important and the 
ground state has two-fold degeneracy. A possible ground 
state can be any superposition of all spin up state |_F|) 
and all spin down state \Fi) where 
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One possible ground state is the GHZ state of N spins 
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both and I GHZ) 
spin entanglement, i.e. 



have the same value of the two- 
zero concurrence, although \P), 
|P|), and \Fi) are separable states but |GHZ)jv is en- 
tangled. The behavior of multi-particle entanglement at 
QPTs is an open problem because a good entanglement 
measure for more than two spins is still under develop- 
ment 0, nil, [3 d . The direct simulation of QPTs via 
AQC might give a clue to this problem. 

For the simulation of QPTs of the XY model via AQC, 
let us decompose the Hamiltonian Hxy into two parts: 
the initial Hamiltonian and the problem Hamiltonian 
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FIG. 1: (a) Phase diagram of the XY chain with the XY 
critical line at A = 1 and the XX critical line at 7 = 0]. 
It is symmetric about 7 and A axes. Path (i) starts from 
the paramagnetic phase and end at the ferromagnetic phase. 
The return path is denoted by (ii). (b) Energy gap Aqi/N 
between the two lowest eigenvalues as a function of 7 and A. 



In the thermodynamic limit, N ^ 00, the ground could 

be either |F|) or |F|) due to the spontaneous symme- jj^^ 

try breaking. However, a computer resource is finite no 

matter whether it is classical or quantum. So we focus 

on the adiabatic quantum simulation of QPTs and the 

generation of entangled states with a finite system. 

The spin 1 /2 XY chain is exactly solvable in the sense 
that its energy spectrum, ground state, and phase di- 
agram can obtained via the mapping of spin operators 
to fermion operators via the Jordan- Wigner transforma- 
tion d, 1^. Recently, this system has attracted much 
attention in the study of the relation between entangle- 
ment and QPTs. Osborne and Nielsen ^ studied the 
entropy of a single spin and the two-spin entanglement, 
where |P|) was taken as a ground state in the ferromag- 
netic phase. Osterloh et al. [31 investigated the two spin 
entanglement of the spin 1/2 XY chain and showed that 
the concurrence as a two-spin entanglement measure ex- 
hibits the characteristic features of QPTs. Note that 
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Among various ways of connecting Hq and Hp, two inter- 
polation schemes, linear and square ones, are considered 
to see whether a proper interpolation could reduce the 
run time of AQC. The linear interpolation is given by 



HxYis,-f) = {l-s)Ho + sHp{j), 
and the square interpolation reads 



(7) 



HxYi-s,7) ={l-s^)Ho + s{2-s)Hpij), (8) 

By comparing of Eqs. (O, and ([H]), one obtains the 
time-dependence of the magnetic field A(s) = (1— s)/s for 
the linear interpolation, and A(s) — (1 — s^)/(2s — s^) for 
the square interpolation. The path of adiabatic evolution 
from the paramagnetic phase to the ferromagnetic phase 
is depicted in Fig. [T] (a). At time s — corresponding 
to the limit A — > cxd, the initial state is given by the 
paramagnetic state, Eq. At s = 1/2, i.e., A = 1, 

the system arrives at the XY critical line. The adiabatic 
evolution ends at s = 1, that is, A = 0. We examine 
which of two states, Eqs. (|3|) and ([5]) is the true final 
state by calculating the fidelity between the GHZ state 
and an evolved state as function of s and 7. Note that 
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Wei et al. [14[] used the maximum fidelity between a state 
and an untangled state as a global entanglement measure 
in the study of the multi-particle entanglement of a XY 
chain. To check the reversibility of AQC, we investigate 
the reverse path from the ferromagnetic phase, starting 
with 11]) or ([5]) , to the paramagnetic phase by exchanging 
and Hp. 

For the numerical simulation of the AQC, we de- 
velop the program which solves the Schrodinger equa- 
tion and diagonalizes the Hamiltonian directly without 
the Jordan- Wigner transformation. We simulate the spin 
1/2 XY chain with TV = 12, and a two-dimensional Ising 
model of size 3 x 3 on a personal computer. 

For the one-dimensional XY model, the energy spec- 
trum is obtained as function of A and the anisotropy pa- 
rameter 7 as shown in Fig. [T] (b). As the free energy 
determines classical phase transitions, the energy gap be- 
tween the ground state and the first exited state plays a 
key role in QPTs. As depicted in Fig. [T](b), the gap Aqi 
between the two lowest eigenvalues of the XY Hamilto- 
nian clearly vanishes at the critical line, i.e., at A = 1. 
On this line the two lowest eigenvalues merge together 
and the ground state becomes degenerate, even though 
the system size is finite. 
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FIG. 2: (a) Aoi and A12 as a function of s at 7 = 1 for 
iV = 8, 10, 12. (b) Aoi and A12 as a function of 7 at A = 0.1 
for N = 8, 10, 12. The linear interpolation is used. 

Fig. [2] (a) plots Aoi and A12, the gap between the sec- 
ond and third lowest eigenvalues, as a function of s at 
7 = 1 for iV = 8, 10, 12. Aqi is nearly independent of 
N. However, A12 is more dependent on TV near the crit- 
ical points Sc = 1/2, i.e., Ac = 1. The energy gap A of 
the spin 1/2 XY chain with infinite size is known to be 
A ~ |A — Ac| Aqi and A 12 can be regarded as the en- 
ergy gaps between the ground state and the first excited 
state in the paramagnetic phase and in the ferromagnetic 
phase, respectively. That is A(A) = Aqi (A) for A > 1 and 



A(A) = Ai2(A) for A < 1. So the critical slowing down 
at the critical point [lB| is due to A12 not due to Agi as 
iV — > 00. Let us examine the universality of Aqi, that is, 
independent of 7. Due to the finite size effect, the region 
of the universality defined by Aqi — does not extend to 
the XX line, i.e., < 7 < 1 for the infinite lattices. The 
region satisfying the universality grows with N as shown 
in Fig. [11(b). 



(a) (b) 




FIG. 3: Energy levels per spin and the fidelity as a functions 
of s for the XY model with N = 12 and 7 = 0.75 (a) for the 
linear interpolation, and (b) for the square interpolation. The 
run time T = 20 is taken. The energy spectrum ranges from 
-1 to 1 and is symmetric about x axis. 

Fig. [3] shows the difference between the linear interpo- 
lation and the square interpolation. The linear interpo- 
lation needs more run time T, than the square one. The 
main reason is that the gap A12 for the square interpola- 
tion is larger at the critical region than that for the linear 
one, so the probability for the transition to the excited 
states is reduced. We find that the square interpolation 
could reduce the run time T of AQC in the case of the 
XY chain. 




FIG. 4: Fidelity -Fghz between the GHZ state and the evolved 
state as a function of A(s) and y. N — 10, T = 20, and the 
square interpolation are used. 

Fig. 0] shows the fidelity between an evolved state and 
the GHZ state, Fghz(s) = |(*(s)|GHZ)|2, as a function 
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of s and 7. The paramagnetic state \P) in Eq. ([3]) at s = 
(A ^ 1) adiabatically evolves to the GHZ state |GHZ)jv 
in Eq. ([5|) at s = 1 (A = 0). Especially, at 7 = 1, this 
result is consistent with Dorncr et a/.'s one [l^. They 
showed that for the Ising chain the paramagnetic state is 
transformed to the GHZ state by slowly decreasing the 
magnetic field A. Although due to the finite size effect, 
the fidelity decreases near the XX critical line, it is almost 
independent of of 7. This is one of the characteristic 
features of QPTs, called the universality. As mentioned 
before, the region of the universality is dependent on the 
number of spins N. 

Let us discuss the reversibility of the QPT. The para- 
magnetic state is adiabatically transformed to the GHZ 
state in the ferromagnetic phase by decreasing A. Does 
the GHZ state evolve adiabatically to the paramagnetic 
state ^ even though there is the energy level splitting 
at the XY critical line? By exchanging Hq and Hp, 
the reverse evolution, the return path (ii) in Fig. [T] (a), 
can be implemented. To examine the reversibility of 
AQC, let us consider Hxy{s,j) — f{s)Ho + g{s)Hp{^) 
where /(s) = 4(s - 1/2)2 ^^^^ g^^^ ^ _4g(g _ 

s = (2 — \/2)/4, the paramagnetic to ferromagnetic tran- 
sition happens. The ferromagnetic to paramagnetic tran- 
sition takes place at s = (2 + V2) /4. We find that in spite 
of the energy level merging and splitting during the jour- 
ney, the paramagnetic state is adiabatically transformed 
to the GHZ state and vice versa as shown in Fig. [5](b). 
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FIG. 5: (a) Two functions /(s) and g(^s) as a function s for the 
round trip, (c) Fghz and Fp as a function of s with 7 = 0.8, 
= 10, r = 200. 



Let us discuss the effect of finite system size on the 
AQC simulation of QPTs and the generation of entangled 
states. Although the system considered here is finite, it 
could show important features of QPTs. First, QPTs 
are characterized by the energy level crossing or avoided 
crossing . Fig. [T] (b) shows energy level crossings at 
A = 1. Second, a XY spin chain exhibits the universal- 
ity, which state that the physics is independent of 7 for 
< 7 < 1. The energy gap in Figs. [T] (b) and the fidelity 
in Fig. m are independent of 7, although the universal- 
ity is not perfect due to the finite size effects. Third, 
Fig. [3] shows the abrupt change in the ground state at 
QPT. One unsolved issue is that the entanglement of the 
ground state in the ferromagnetic phase. If either \F^) or 



\F\) is the ground state in the ferromagnetic phase, there 
is no entanglement. On the other hand, the ground state 
|GHZ)Ar is a multi-particle entangled state. In this case, 
the single-site entropy varies from at the paramagnetic 
phase to 1 in the ferromagnetic phase. This is in con- 
trast with the result of ;4j . A related question is whether 
the entanglement of the ground state in ferromagnetic or 
paramagnetic phase is constant or dependent on A. In 
other words, does entanglement change abruptly only at 
the critical point? 



IV. ENTANGLEMENT OF 
TWO-DIMENSIONAL ISING MODEL 

Let us turn to the two-dimensional Ising model to pro- 
duce the two-dimensional GHZ state 



|GHZ)2D = 




(9) 



where is the spin-up state at the lattice site i and 
2- As illustrated in Fig. [5] (a), the 3x3 two-dimensional 
lattice is considered. The open boundary condition is as- 
sumed. The two-dimensional Ising model can be mapped 
to the one-dimensional Ising model with long-range inter- 
actions as depicted in Fig. [S](a). Fig.[n](b) and (c) show 
several lowest energy levels, and the fidelity between an 
evolved state and the two-dimensional GHZ state defined 
by Eq. ([H]) as a function of s. Like one-dimensional spin 
1/2 XY model, the two dimensional GHZ state can be 
generated via AQC. 

(a) 





(c) 
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FIG. 6: Two dimensional square lattice of size 3x3 and 
its mapping to the one-dimensional lattice with long-range 
interaction, (b) Energy levels as a function of s with the 
hnear interpolation, (c) Fidelity between the evolved state 
and the two-dimensional GHZ state as a function of s. 



V. SUMMARY 

We have considered the one-dimensional XY model 
and the two-dimensional Ising model and simulated the 
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QPTs and the generation of an entangled state via AQC. 
Although the system size is finite, our results show the 
characteristic features of QPTs. It has been demon- 
strated that the paramagnetic state evolves adiabatically 
to the GHZ state in the ferromagnetic phase and vice 
versa. The generation of entangled states via AQC is 
simple in the sense that only the external magnetic field 
is turned off or on slowly. It doesn't require the control 
of the exact qubit-qubit coupling, i.e., CNOT gate. We 
have shown that a square interpolation scheme is better 
in reducing the run time than linear one. 

One open issue in AQC is that the minimum energy 
gap, which determine the run time T, should be known 
before running. Also the run time should be smaller than 
the decoherence time but at the same time large enough 



to avoid the unwanted transition. We are studying the 
effect of decoherence on AQC by solving the Lindblad 
master equation and the generation of W-type en- 
tangled states or cluster states [l^ via AQC. Also, it is 
interesting to study how to simulate a quantum system 
in the thermodynamic limit with finite quantum compu- 
tational resources. 
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